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Abstract 

In a model of large distance modified gravity we compare the nonperturba- 



tive Schwarzschild solution of hep-th/0407049 to approximate solutions ob- 



tained previously. In the regions where there is a good qualitative agreement 
between the two, the nonperturbative solution yields effects that could have 
observational significance. These effects reduce, by a factor of a few, the pre- 
dictions for the additional precession of the orbits in the Solar system, still 
rendering them in an observationally interesting range. The very same effects 
lead to a mild anomalous scaling of the additional scale-invariant precession 
rate found by Lue and Starkman. 



1 Introduction 



The DGP model of large distance modified gravity [T| has one adjustable parameter 
- the distance scale r c . Distributions of matter and radiation which are homoge- 
neous and isotropic at scales > r c exhibit in this model the following properties: for 
distance/time scales <C r c the solutions approximate General Relativity (GR) to a 
high accuracy, while for scales > r c they dramatically differ [H El El El ■ Postulat- 
ing that r" 1 ~ H Q ~ 10~ 42 GeV the deviations from GR could lead to interesting 
observational consequences in late-time cosmology, see, e.g., [HIE], |5]-[TT]. 

On the other hand, sources of matter and radiation with typical inhomogeneity 
scale less than r c have somewhat different properties. These are easier to discuss for 
a Schwarzschild source - a spherically-symmetric distribution of matter of the mass 
M and radius r , such that r M < r < r c [r M = 2G^M is the Schwarzschild radius 
and GV the Newton constant). For such a source a new scale, that is a combinations 
of r c and r M , emerges (the so-called Vainshtein scale 1 ) jl]: 

r^(r M r c 2 ) 1/3 . (1) 

Above this scale gravity of a compact object deviates substantially from the GR 
result. Note that r* is huge for typical astrophysical objects. An isolated star of a 
solar mass would have r* ~ 100 pc. However, if we draw a sphere of a 100 pc radius 
with the Sun in its center there will be many other starts enclosed by that sphere. 
The matter enclosed by this sphere would have even larger r*. We could draw a 
bigger sphere, but it will enclose more matter which would yield yet larger r* and so 
on. An isolated object which could be separated from a neighboring one by a distance 
larger than its own r* is a cluster of galaxies. For typical clusters, r* ~ (few Mpc) 
is just somewhat larger than their size and is smaller than their average separation. 
The above arguments suggest that interactions of isolated clusters will be different in 
the DGP model. On the other hand, at scales beneath a few Mpc or so, there will be 
agreement with the GR results with potentially interesting small deviations. Below 
we discuss these issues in detail on an example of a single isolated Schwarzschild 
source. There exist in the literature two different solutions for the Schwarzschild 
problem in the DGP model. The first one is based on approximate expansions in 
the r<Cr t and r ^> r* regions [UIHII3] ( see a l so [111 USD- We call this set of results 
the perturbative Schwarzschild (PS) solution. The second one [TH] is a solution 
that interpolates smoothly from r « r, to r > r c > r», and is non-analytic 
in the either parameters used to obtain the PS solution. We call this the non- 
perturbative Schwarzschild (NPS) solution. It is important to understand which of 
these two solutions, if any, is physically viable. Since neither of the two can be solved 
completely without numerical simulations, a first step to discriminate between them 
would be to look closely at the theoretical differences, as well as predictions that 
could by tested observationally. This is the goal of the present note. 

*A similar, but not exactly the same scale was discovered by Vainshtein in massive gravity |12) . 
hence the name. 
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2 Qualitative discussions 



We will study separately two regimes, r<Cr t and r ^> r*. 

(I) r<r t . In this regime the standard Gn expansion breaks down jl]. How 
could one proceed? One way is to perform an expansion in powers of m c = r~ l 
[3]. This expansion breaks down above r ~ r* but is well suited for the r « r, 
domain (Kaloper J7| recently used a different expansion. His proposal could prove 
to be useful for a broad class of problems). A Schwarzschild metric in the small 
m c expansion was calculated by Gruzinov [13] (see also [El)- ft i s instructive to 
compare the result of ^3] with the NPS solution of [To] . 

Let us start with the Newton potential <f>(r). The expansion of the exact result 
of [in] for r <C r* leads: 



where = 3/2- 2( v / 3 - 1) ~ 0.04, and a is a number to be discussed in detail 
below. The above result, but with (3 = 0, is what was first obtained in a small m c 
expansion jT3j. The NPS solution of ^E] gives (3 ~ 0.04, it depends on irrational 
powers of m c jTH] , and it differs by that from the small m c expansion results. 

Is the above difference important? As was demonstrated in Refs. ^3] and 
[IS], the modification of the Newton potential in (J2J, although tiny, could lead to 
a measurable precession of orbits in the solar system (see, Refs. ^H] for further 
studies). The above works used the potentials obtained in the small m c expansion, 
e.g., used (J2J with {3 = 0. Although {3 is tiny, the ratio (r*/r) is typically huge in 
the cases of interest, therefore, taking into account the effects of a nonzero (3 could 
lead to appreciable differences in the predictions of the PS and NPS solutions. We 
will study this issue in the next section. 

Consequences of the modified potential (j2J) could be understood as well in terms 
of invariant curvatures. The Schwarzschild solution in GR has zero scalar curvature. 
In contrast with this, the solution (J2J) generates a nonzero Ricci scalar that extends 
to r ~ r, in the NPS solution (see, [IH] and discussions below). This can be seen by 
looking at the trace equation in the DGP model: 



where R is the 4D Ricci scalar, K is a trace of an extrinsic curvature and T is a 
trace of the stress-tensor times 87tGn (for the ADM formalism in the DGP model 
see, e.g., [213 EI])- This has to be compared with the trace equation in GR: R = T. 
The second term on the LHS of (jSJ) is not zero outside the source and, therefore, 
gives rise to nonzero R. This curvature, although tiny, extends to enormous scales 
of the order of r ~ r* [16 . The sign of the curvature depends on a choice of the 
boundary conditions in the bulk, since the latter determines the sign of K. There 
are two choices for this. The so-called conventional branch corresponds to a negative 





R — 3m c K = T , 



(3) 
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(AdS like) curvature produced by the Schwarzschild source, while the selfaccelerated 
branch [2] corresponds to a positive (dS like) R. This is reflected in the sign of the 
coefficient a in which takes a positive value on the conventional branch and 
becomes negative on the selfaccelerated branch: a ~ ±0.84. Therefore, there is 
an additional tiny attraction toward the source on the conventional branch and a 
repulsion of the same magnitude on the selfaccelerated branch. This change of sign 
was first found by Lue and Starkman ^3] in the context of the PS solution. 

(II) r ^> r*. In this regime the small m c expansion breaks down. However, the 
conventional Gn expansion can be readily used [HE]. The results are [Tj: 

(A) For r ^> r* DGP gravity is a tensor-scalar theory, where the extra scalar 
couples to matter with the gravitational strength: the vDVZ phenomenon fZ*2\ \'23\. 

(B) The Newton potential scales as 1/r for r» C r C r c which smoothly transi- 
tions into the 1/r 2 potential at r ^> r c . 

These properties of the PS solution were reconfirmed in detailed studies of Refs. 
[13 UH 113 123 123- Could the PS solution interpolate from r « r, to r » r c ? The 
above question is related to the following one: what is a gravitational mass that is 
felt by an object separated from the source at a distance r ^> r*? The PS solution 
implies that this is just the bare mass M of the original source. On the other hand, 
one may expect that the curvature created by the source in the domain r <C r* 
would also contribute to this effective mass (the ADM mass) [16j. If so, unless there 
is a hidden nontrivial cancellation, a putative observer at r > r t would measure 
an effective mass different from M. The above property is captured by the NPS 
solution of Ref. It has the following features: 

(A') For r 3> r* it is a solution of a tensor-scalar gravity (as in (A) above); 

(B') The Newton potential scales as 1/r 2 for r ^> r* (different from (B)). 

An attractive feature of the NPS solution is that it smoothly interpolates from 
r < r„ tor >> r c . However, a somewhat unusual fact is that it does not recover the 
results of the Gn expansion. This will be discussed in the reminder of this section 
(readers who are not interested in these somewhat technical issues could directly go 
to the next section without loss of clarity). 

Why is that, that the NPS solution (THj does not agree with the results of the 
perturbative Gn expansion, even in the regime r ^> r*, where the latter approxi- 
mation is internally self-consistent? There could be a few different reasons for this. 
Formally, one is solving nonlinear partial differential equations and these can have 
different solutions even with the same boundary conditions. In our two cases, how- 
ever, the boundary conditions are somewhat different: the PS solution is supposed 
to describe the same mass M at short and large distances, while the NPS solution 
matches M at the short scales but asymptotes to a screened mass at the large scales 
2 . Then either the PS and NPS solutions belong to different sectors and are both 
stable, or at least one of them should be unstable. In the former case, one should 
distinguish between them observationally, while in the latter case a relevant point 

2 The boundary conditions at the brane are also different, see a footnote on page 6. 
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would be that the ADM mass of the NPS solution is smaller [TJJ]. In a very qual- 
itative way, this can be understood as follows. A deviation from the conventional 
metric at r < r, sclaes as m cy frMr (we ignore small f3 here.) This can give rise 
to a scaling of the scalar curvature wicy/rur'^/ 2 . The curvature extends roughly to 

r ~ r*, and the integrated curvature scales as m cv /rj\7^*^ 2 ~ r M . Then, the "effec- 
tive mass" ' due to this curvature can be estimated as r^/Mpj ~ M, which is of the 
order of the mass itself. 

On the other hand, it may well be that there is a certain "discontinuity" between 
the linearized and full non-linear versions of the DGP model in 5D. This could result 
from a different number of constraints one has to satisfy depending on whether 
solutions are looked for in the linearized approximation or in the full non-linear 
theory. For instance, one of the bulk equations can be combined with the junction 
condition in 4D to yield: 

Zm\R = R 2 - 3Rl„ . (4) 

On a flat background both terms on the RHS of contain at least quadratic 
terms in the fields. Therefore, according to R has to be zero in the linearized 
approximation. The latter condition happens to be a consequence of the other 
linearized equations of the theory as well; therefore, (j3J) is trivially satisfied as long 
as those other equations are fulfilled. This changes at the nonlinear- level: Eq. 
(0J) becomes an additional constraint that one has to satisfy on top of the other 
equations. Because of this: (i) The solutions of the linearized theory may not be 
supported by the nonlinear equations (the phenomenon known as "linearization 
instability" in gravity), (ii) New non-perturbative solutions that do not exist in 
the linearized theory may emerge. One way to decide on the point (i), is to study 
solutions for other sources and see whether a similar phenomenon takes place. The 
NPS solution of 16 is an explicit example of the point (ii). 

3 Explicit solution 

We consider the action of the DGP model pQ: 

S = Ml J d 5 x^/=g~R + M 2 P j d 4 x^-^R . (5) 

Here, the (4 + 1) coordinates are x M = (x^, y), fi = 0, . . . , 3 and g and R are the 
determinant and curvature of the 5 dimensional metric QmNi while g and R are 
the determinant and curvatures of the 4 dimensional metric g^ = g /Jtl/ (x' J ',y = 0). 
The Gibbons- Hawking surface term that guaranties correct equations of motion 
is implied in the action (0). Mp denotes the 4D Planck mass and is fixed by the 
Newton constant. On the other hand, the scale M* is traded for the parameter 
r c = M 2 /2Ml discussed in the previous section. 
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The NPS solution studied in [THj is found by considering a static metric with 
spherical symmetry on the brane and with Z 2 symmetric line element: 



ds 2 



-e~ x dt 2 + e x dr 2 + r 2 dVL 2 + 7 drdy + e a dy 2 



(6) 



where A, 7, a are functions of r = ^/x^g^ u x u and y. The Z 2 symmetry across the 
brane implies that 7 is an odd function of y while the rest are even. The brane is 
chosen to be straight in the above coordinate system 3 . 

The exact solution for y —> 0+ is given implicitly as follows: 



1 



P(r) 



where P is obtained from 



3 r 

P(r) = —-^m 2 c J dr r 2 U(r) , 



(7) 



(8) 



in which U can have two different behaviors corresponding to the solution of the 
following two equations (giving rise to a conventional and selfaccelerated branch 
respectively) : 



(k 2 rf 



;i+3c/+/) 



U 2 (3 + 3U + y/3f) 2 ^(-5 -SU + f) 
(-5 -3U + /)(-3 -3U- V3f) 2V * 



(U + 2) 



3U + f) 



where / = yl + QU + 3U 2 and k is an integration constant. 

Note that in this parametrization the gravitational potential 
approximation is easily obtained, namely 

P(r) 



2r 

The off-diagonal metric component, 7, is determined from 

2r 2 P r _ ( r4 7e~ A ) r 



(9) 
(10) 

in weak field 
(11) 



Pr 



( r7e A) 



(12) 



and the yy component from 



2r 2 P r . 



+ e" V 



(13) 



3 One could transform ^ to the coordinate system where the metric is diagonal ds 2 - 
-A(r, z)dt 2 + B(r, z)dp 2 + C(r, z)dfl 2 + dz 2 , and A ^ B. In this system our brane will be bent. 
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The profile X y for y — > 0+ can be computed as well: 

A y = e" A 7r . (14) 

The two integration constants, k and the one produced in the integration (|8j). are 
determined by imposing appropriate boundary conditions at the source (namely, 
P(r —>■ 0+) — > tm) and at large distances, (namely, A ~ ^|//^ 2 in the conventional 
branch or A ~ m 2 r 2 + rf^/r 2 in the selfaccelerated branch and no 1/r term). 



3.1 Conventional branch 

The conventional branch is obtained from the solution of Q- As shown in ^E] 
the boundary conditions (-P(O) = r^, P(+oo) = 0) determine the exact relation 
between k\ and r*, namely 

2(nk) 3 = c , (15) 

where c is the following integral: 



[7 2 (3 + 3C/ + v / 3/) 2V5 (-5 - 317 + /) 



dC/ « 0.43 . (16) 



The solution has the following asymptotic behavior. At large distances, r ^> r* 
([/ — > + ), we obtain 



P(r) -^+..., (17) 



where, 



3V2 



r ". = ipWlp¥" a56rMr " (18) 



2 



while at short distances, r C r t ({/ ^ +00), we get 



P (r)_r M ai ^(rA^\..., (19) 



iy iy \ 7" 

where 

As we see, a short distance observer at r^f <C r <C r* would measure the gravitational 
mass M with a small corrections to Newton's potential, while the large distance 
observer at r > r t would measure an effective gravitational mass ~ M(jm /t c ) 1 ^ 
[TE] . The latter includes the effects of the 4D curvature. 
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3.2 Selfaccelerated branch 



The solution on the selfaccelerated branch is obtained from The relation 

between k\ and r* is obtained, as in the conventional case, by imposing boundary 
conditions (-P(O) = tm, P(r) — m 2 r 3 — ► for large r). This gives 



£/ 2 (3 + 3C/ + \/3/) 2V 3(-5 — 3Z7 + /) 



dC/ 



= 6 3v ^ /4 c « 4.41 . (21) 

The second line in (I21J1 . that is generated by a change of variables in the integral 
(U = —U — 2) while using (fTKj) . also gives a relation between fcj. and &2, 

fc 2 = 6^/^! . (22) 

The solution has the following asymptotic behavior. At large distances, r ^> r* 
(U — > — 2 _ ), we derive 

P(r) --^ M2 +m c 2 r 2 + ... , (23) 



where, 
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f M 2 = ~ 7^ if « °- 45 V, , (24) 



(3 - v^) 2 ^ ^2 
while at short distances, r <r, (Z7 — > — oo), we get 



- a 2 m c r^ ( "T ) + • • • , (25) 



where a 2 = — cki ~ —0.84 is, in absolute value, the same constant appearing in the 
conventional branch short distance expansion ([19)1 . Note, however, that the sign of 
the correction to the AD behavior is opposite in the two branches. 

At intermediate distances, r„ < r < r c , the potential contains a 5.D gravitational 
term that is repulsive, f 2 M /r 2 . This looks like a 5D negative mass. However, this 
is not an asymptotic value of the mass since one can only cover the solution in 
the above coordinate system till r ~ r c where the dS like horizon is encountered. 
Moreover, in the intermediate regime r t <r < r c , the de Sitter term m 2 r 2 in the 
potential always dominates over the f 2 M /r 2 term suggesting that the effects due to 
the Schwarzschild source are strongly suppressed. 



3.3 Perihelion precession 

The deviation from 4D gravity (J2J gives rise to the additional perihelion precession of 
circular orbits [TJJ [TH] (see also |TH] for comprehensive studies of these and related 
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issues). In a simplest approximation this effect is quantified by a fraction of the 
deviation of the potential from its Newtonian form 



e=-f. 26 

This can be used to evaluate an additional perihelion precession of orbits in the 
Solar system [HUTS] 4 . As we discussed in Section 1, the e ratio is somewhat different 
for the non-perturbative solution (NPS solution) as compared to the approximate 
solution (the PS solution ) used in Refs. [JJJ HHj. We can easily calculate this 
difference: 



(-NPS \ a \ 



/ r \ P / r \ 0.04 

(") -0.59 (-) . (27) 



e PS V2 
The perihelion precession per orbit is 



3nr M 37r|d /r\ 3 / 2 / r \ 004 
A^ = 2tt + — =F — - — - ( — ) - . (2f 



The second term on the RHS is the Einstein precession, and the last term arises 
due to modification of gravity. For the PS solution this was first calculated in Refs. 
f2JE|; the solution (J28|) is written for the NPS solution and is somewhat different. 

For the Earth-Moon system r = 3.84 x 10 10 cm and rf arth ~ 6.59 x 10 12 cm; as 
a result the ratio in ()27|) is approximately 0.48. Therefore, the predictions of the 
NPS solution for the additional perihelion precession of the Moon is a factor of two 
smaller than the predictions of the approximate solution. The result of (|2*H|) for the 
additional precession (the last term on the RHS) is =f0.7 x 1CT 12 (the plus sign for 
the selfaccelerated branch). This is below the current accuracy of 2.4 x 1CT 11 [2*7] . 
but could potentially be probed in the near future [21] 5 - 

A similar calculations can be performed for the anomalous Martian precession 
|14[ I18j. For the Sun-Mars system we get: 

/ \ 0.04 

— ~ 0.59 1 TSun : Mars ) - 0.30, (29) 
eps V r * un J 

where we used rs un -Mars — 2.28 x 10 13 cm and rf un = 4.9 x 10 20 cm. Therefore, we 
see that the suppression in the NPS result for the precession of the Martian orbit 
is stronger. The additional precession of the Mars orbit is ~ =pl.3 x 10 -11 , which 
should be contrasted with a potential accuracy of the Pathfinder mission ~ 9 x 10~ n . 

4 Note that in the leading order of the relativistic expansion the answer is given by the correction 
to the Newtonian potential, while the correction to the rr component of the metric is not important. 

5 An interesting possibility that similar effects could leed to seemingly observable increase of the 
Astronomical Unit was recently discussed in |29| . 
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Last but not least, Lue and Starkman (LS) [Hj, found that the PS solution gives 
rise to a correction to the precession rate (additional precession per unit time), 




(30) 



that is universal, i.e., is independent of the source. The NPS solution, predicts a 
weak anomalous violation of the universal Lue-Starkman scaling due to the RHS of 
(J27J) . The results is 



This rate depends mildly on the source mass and a separation from it. The rate is 
a slowly increasing function or r, as opposed to the rate due to the second term on 



the RHS of (J28|) . which is decreasing with growing r as r Einstein = y^r^/Sr 5 . 



In this note we compared the PS [TJ EJ HH1 EH HHj and NPS |TE] solutions in the 
DGP model. We emphasized different, but interesting predictions that these two 
solutions make in the observationally accessible domain of r <C r*. These predictions 
are testable. 

As we have also mentioned, there will be important differences in the predictions 
at r > r,. These need further detailed studies, especially in the context of the 
structure formation. We would expect that both the linear as well as non-linear 
regimes of the structure formation will be affected. If the NPS solution is the 
right one, then even at very large scales nonperturbative techniques should be used. 
Moreover, the nonlinear regime of the structure formation could be sensitive to, and 
be able to discriminate between, the PS and NPS solutions. 

The same issue of nonlinear interactions arises in the context of strong coupling 
behavior in the 5D DGP model |U E01 ED E21 E3] • This is related to the problem of 
the UV completion of the quantum theory [3U1 EI] for which seemingly two different 
proposals were put forward in Refs. [32] and [33] ■ It would be interesting to pursue 
these studies further. The string theory realizations of brane induced gravity of Refs. 
|3H ESI EE] can be taken as a guideline. It would also be interesting to understand 
the NPS solution in terms of the approach of Refs. [3U1 E3J. 

We have not touched upon the issue whether the small fluctuations on the self- 
accelerated branch contain a negative norm state [3UJE3], or not ( see also |57|). 
and when these fluctuations are relevant. Additional investigations on this issue are 
being conducted. 

It would also be interesting to look at the Schwarzschild solutions in models of 
large distance modified gravity where nonlinear interactions do not exhibit the strong 
coupling behavior. This is the case [3S] in a certain models of brane induced gravity 




(31) 




4 Outlook 
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in more than five dimensions EE] > as well as in the "dielectric regularization" of 
the 5D DGP model j3U]- Finally we would also point out the constrained approach 
to the 5D DGP model [¥H W]\ in which case strong interactions also seem to be 
absent. All the above deserves further detailed investigations. 
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